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The curved space-time around current loops and solenoids carrying arbitrarily large steady electric 
currents is obtained from the numerical resolution of the coupled Einstein-Maxwell equations in 
cylindrical symmetry. The artificial gravitational field associated to the generation of a magnetic 
held produces gravitational redshift of photons and deviation of light. Null geodesics in the curved 
space-time of current loops and solenoids are also presented. We hnally propose an experimental 
setup, achievable with current technology of superconducting coils, that produces a phase shift of 
light of the same order of magnitude than astrophysical signals in ground-based gravitational wave 
observatories. 
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I. INTRODUCTION 

Somehow, studying gravity is a contemplative activity: 
physicists restrict themselves to the study of natural, 
pre-existing, sources of gravitation. Generating artificial 
gravitational fields, that could be switched on or off at 
will, is a question captured or left to science-hction. 

However, the equivalence principle, at the very heart 
of Einstein’s general relativity, states that all types of 
energy produce and undergo gravitation in the same 
way. The most widespread source of gravitation is the 
inertial mass, which produces permanent gravitational 
fields. At the opposite, electromagnetic fields could be 
used to generate artificial, or human-made, gravitational 
fields, that could be switched on or off at will, depending 
whether their electromagnetic progenitors are present or 
not. 

The equivalence principle actually implies that one 
also generates gravitational fields when generating 
electromagnetic fields. However, since the gravitational 
strength is extremely small compared to the one of the 
electromagnetic force pT]. large electromagnetic fields 
will only produce tiny space-time deformations. Yet, 
electromagnetic fields do curve space-time. Therefore, 
general relativity predicts that light and more generally 
electrically neutral massive particles are deflected by 
electromagnetic fields, although they do not feel the 
classical Lorentz force. This effect does not require new 
exotic physics and it might serve in the future to build 
new tests of the equivalence principle in the laboratory. 
In experimental gravity, the permanent gravitational 
fields involved cannot be withdrawn completely. At 
the opposite, the gravitational fields generated by 
electromagnetic fields can be switched off: their ex¬ 
perimental search can therefore be done by comparing 
measurements made in presence and in absence of 
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electromagnetic fields. However, due to the weakness 
of the gravitational interaction, even the strongest 
magnetic fields humans can currently generate will only 
produce tiny space-time deformations. Detecting them 
would constitute a true experimental challenge which 
we glimpse at in this paper. Such a detection would 
nevertheless open the way to new laboratory tests of the 
equivalence principle. 

The so-called Einstein-Maxwell (EM) equations re¬ 
group the classical field theories of general relativity 
and electromagnetism in a covariant way, although 
without truly unifying them. In some sense, the idea of 
gravitational field generation from a magnetic field can 
be attributed to Levi-Civita whose early analytical work 
[1] describes the curvature of space-time completely filled 
by a uniform magnetic field. Subsequent works have 
established analytical solutions for space-time around 
an infinitely long straight wire carrying steady current 
[31]. The main problem of these analytical solutions 
is that the associated metric is not asymptotically fiat, 
due to the infinitely large current distribution, which 
makes these analytical solutions of poor interest for 
practical applications. Overcoming this problem requires 
considering current distributions of finite extent such 
as current loops and solenoids. Asymptotic space-time 
around current loops carrying steady current has been 
studied in |4]. An attempt to derive the full solution of 
EM equations around the current loop was realized a 
bit later [5], however this attempt lead to an unphysical 
solution due to an oversimplifying assumption. The case 
of an infinitely long solenoid was considered in [6] but 
only for weak perturbations of the metric in linearized 
general relativity. Therefore, the solutions of the full 
non-linear EM equations sourced by the steady currents 
carried by loops and finite solenoids remained so far 
unexplored until now. 

In this paper, we present two important results: 
(1) how space-time is curved around current loops and 
solenoids carrying arbitrarily large electric currents 
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and (2) how the consequent deviation of light could be 
detected. The structure of this paper is as follows. In 
section II and III, we develop the numerical resolution of 
coupled EM equations in cylindrical symmetry. We then 
present the trajectories of light in curved space-times 
around current loops and solenoids in section IV. In 
section V, we finally propose an experimental set-up, 
based on a modification of interferometers used for the 
search for gravitational waves and current technology 
of superconducting electro-magnets, that produces a 
detectable artificial space-time curvature. We conclude 
by emphasizing the importance of the effect presented 
here: the deflexion of light in the curved space-time of 
an electromagnet opens the way to new experimental 
tests of Einstein’s equivalence principle. 


II. EINSTEIN-MAXWELL EQUATIONS FOR 

THE CURRENT LOOP AND THE SOLENOID 

A. Field equations 

The EM system models the interaction of gravitation 
and electromagnetism by their juxtaposition in the fol¬ 
lowing coupled tensorial field equations (in S.I. units [18]): 

_ (em) 

~ fll' \^J 

= /ioJ" (2) 

where = -^ is the 

Maxwell stress-energy tensor, Ay — dyA^ is the 

Earaday tensor of the electromagnetic field, g^y and A^ 
are the metric and the four-vector potential, i.e. the fun¬ 
damental fields describing gravitation and electromag¬ 
netism respectively, R^y is the Ricci tensor and the 
four-current density. 

Space-time is therefore curved by the energy of the 
electromagnetic field as ruled by Einstein equations of 
general relativity Q. In the same time, the electro¬ 
magnetic field propagates in the non-trivial background 
it generated through Einstein’s equations Q, and this 
propagation is described by the covariant Maxwell equa¬ 
tions in curved space-time ([^ . Because we are interested 
in the additionnal gravitational field that is produced by 
the magnetic field, we neglect the mass of the current 
carriers and the electric wires [T9|. 

Since current loops and solenoids possess one axis of 
symmetry, we choose the so-called Weyl gauge laisiiz] 
for the metric field: 

ds^ = + dz'^) - 

( 3 ) 

In this symmetry, the vector potential A^ trivially re¬ 
duces to one non-vanishing magnetic component A^ = 
a{r^z)lr. The advantage of the Weyl gauge @ is that the 
equations of motion directly exhibit usual Laplacian op¬ 
erators on fiat background with cylindrical coordinates. 


Indeed, Eqs. 0 and 0 using Eq.(|^ now read (see also 
0): 


d^X + d^p 


8ttG 

c‘^po 

SnG 


{{draf - {dzaf^ 


rdrpdzp + 


IGttG 
C^po r 


drCidzCi 


( 4 ) 

( 5 ) 

(6) 


- (drttdrp + dzttdzp) - rpoJ (7) 


where + 3“^ the usual Laplacian on 

flat space in cylindrical coordinates and where J is the 
angular component of the current density. In [S], the 
term (drttdrp + dzad^p) in Eq. 0 was arbitrarily set to 
zero to provide an analytical solution. We do not restrict 
ourselves here to this arbitrary reducing condition and 
provide the full solution numerically. 

Eq.g is the Maxwell equation on a curved space- 
time described by cylindrical coordinates. Eor a flat 
Minkowski background p = A = 0, we have that the 
non-relativistic field Unr satisfies 

2 

~^r^nr ~ (8) 

We can therefore decompose the total field a(r, z) into 
the sum of a non-relativistic part Unr (solution of Eq.Q) 
and a relativistic contribution Urei by setting a = Unr + 
Urei- The relativistic contribution Urei will therefore be a 
solution of 

^(r,2:)^rel —drOj^el — {(^r^nr T ^r^rel) T • * • 

(a^Unr + dza^el) d^p} (9) 


This avoids dealing with point-like sources representing 
the current loop and the solenoid in cylindrical coordi¬ 
nates. The source of the field equations now lies in the 
non-relativistic contribution anr(^, z). 

A current loop of radius I corresponds to a current den¬ 
sity located on a infinitely thin ring: J ^ S{z) • 6{r — 1) 
while a solenoid of finite length L and of radius I corre¬ 
sponds to a current density located on an infinitely thin 
sheet located at r = / and 2 : G §] [H]- Analytical 
expressions of the vector potential in both cases can 
be derived from the Biot-Savart law, expressions that of 
course verify Eq.([^. The non-relativistic solution 
for the current loop is given by BM- 

(i») 

where / is the steady current carried by the wire, 

= Frl ((/ -1- r)^ + z‘^) 








and 
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K{k^) = [ (l —sin^((/:?)) dip, 
Jo 


g^ = 


Arl 


(/ + r)2 


/*7r/2 

E{k^) = / (l —sin^((/:?)) 
Jo 


1/2 


dip 


are the complete elliptic integrals of the first and second 
kind respectively. 

For the solenoid of finite length L, the non-relativistic 
solution z) is given by EDI: 






k‘^g^ 


E{k^) 5^ -1„, 2 


k‘^ 




( 11 ) 
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where n is the number of wire loops per unit length, 


{{l + ry+e)’ 


U = z±^ 


and 


n(y,fc^)=/' (l-y sin2((^)) yi - fc^sin^(</?)) 
Jo 


■ 1/2 


dip 


is the complete elliptic integral of the third kind. 

In order to be used as source terms in Eq.(|^, the gra¬ 
dients of anr (!/_£) obtained analytically from the 

formulae Eqs( 1Q][TT ) and the properties of complete el¬ 
liptic functions. 


If we set r = ul, z = vL and anr,rei ^ ^inr,rei/ ik^oE) (anr,rei ^ <i^nr,rei/ (ko'^EL) for the solcuoid), Eqs. ( 4|5|9 ) now 
reduce to the following set of dimensionless equations 


j^2 gp / j2 

V p — f (^u(^nr H“ ^rel)) H“ (^nr H“ ^rel)) 


f 


eP 


P V? 


V A + ^ ^ ( ((9«(anr + Orel))-(Onr + Orel)) 


L2 


2 / /2 
V ttrel “^u^rel — f ^u(^nr H“ (^ve\)duP “1“ H“ ^rel)^vP 


0 — uOupdyp dyP “|“ 2iCj ^n(^nr “t“ ^lrel)^v(^nr H“ ^rel) 


( 12 ) 

(13) 

(14) 

(15) 


where = d ^, E ^du + j 2 dy. In the following, we will solve Eqs. ( 12|l4) numerically and use the last of the Einstein 
equations Eq.(15) for validation check. _ 


The dimensionless magneto-gravitational coupling for 
the current loop and the solenoid are given by 

; cf = (16) 

Hence it is the square of the total current (/ for the loop 
and InL for the solenoid) that sources the gravitational 
field [20] 


or the solenoid {{u,v) (+oo,±oc)), these devices be¬ 

have as magnetic dipoles and the space-time is asymp¬ 
totically flat (see also 0). The magnetic component 
a is then ruled by Eq.Q with dudduP + dyddyp ^ 0 
(i.e., the non-relativistic Maxwell equation). This condi¬ 
tion is achieved if Urei vanishes and S^anr = udu"^ and 
dudnv = —udy^jj with 'll; an harmonic function called the 
scalar magnetic potential H E]. Eor magnetic dipoles, 
the (dimensionless) scalar magnetic potential is given by 


B. Boundary conditions 


Solving the system ( T2p4 ) requires the specification of 
boundary conditions. On the axis of symmetry, r = 0, 
space-time must be smooth such that we have drp\r=o = 
drX\r=o = drd\r=o = 0- Far away from the current loop 


The metric functions p and A are therefore ruled by 
the following equations at large distances ((u,v) 
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(+ 00 , ±oo), p ^ 1 so that ^ 1): 

vv = Ci^ ^ 

v^A = Ci^ {d,^f - {du'iPf^ ■ 


The asymptotic behaviors for the metric fields around 
the current loop and the solenoid are therefore given by 


P 

A 


Cx L‘ 2 

32 ’ 

CxL^ 


1 


-3 


16 P 


2u^ 


2 

“ + p” 


-3 


L ‘^ f 2 ^^2 

■PTl” +■?” 


-3 




-4' 


(18) 


(19) 


III. NUMERICAL RESOLUTION OF FIELD 
EQUATIONS 

A. Numerical Method 


We solve Eqs.( T2][l4 ) numerically by using a combi¬ 
nation of relaxation and spectral methods. We first in¬ 
troduce a sequence of functions X^^\u,v) and 


a['^}(u,v) for the relaxation algorithm such that Eqs.(12 
|14| ) can be approximated by a set of linear inhomoge¬ 
neous elliptic equations: 


V2p("+1) 

V2a(«+1) 


y 2 (n+1) _ 2 „ (n+1) 

®^rel 


51 

52 

53 




P > “rel 
«(") 


( 20 ) 

( 21 ) 

( 22 ) 


where the source terms Si gather all non-linear terms of 
Eqs.( T2][T4 ) but evaluated with the previous state n of the 
relaxation algorithm. The algorithm starts with a state 
n = 0 that corresponds to the non-relativistic solution: 
p^^\u^v) = {u, v) = 0. At each relax- 


ation step, we solve Eqs.([^[^ with a spectral method 
to compute and be¬ 

fore iterating. The algorithm is stopped when the rela¬ 
tive update of the three fields , averaged over the spatial 
domain, reaches some tolerance threshold. 

To solve Eqs.( 2Qp2 ) at each relaxation step, we de¬ 
velop each field fjupv) and each source term Si{u,v) as 
a truncated Eourier series in the 'i;—direction: 


since all the fields are even functions of v due to cylin¬ 
drical symmetry and where 


fk{u) = 4 J ^f{u,v) 


) cos ( 1 dv- 


In Eourier space, Eqs.( 2Qp2 ) now become a set of linear 
inhomogeneous ODEs 


(fpk 1 dpk P (fcTrV , 
du^ ^ u du L'^yv) 

d'^Xk 1 dXk P (\ 

du^ ^ u du L^\v) 

d^gpelfc I dUrelfc P (k'K\^ . 

du^ u du \ V J 


5fc,i(u)(23) 

5fc,2(u)(24) 

Sk,3N(25) 


for k = 1, • • • , A and where we have omitted the relax¬ 
ation index n + 1. Eqs.(^25 ) can be solved as a bound¬ 
ary value problem whose boundary conditions at n = 0 
and u = U ^ 1 are the Eourier transforms of the con¬ 
ditions derived in section II.B. In the present paper, we 
have used the algorithm described in m for the numer¬ 
ical resolution of Eqs.( [23p5 ) for each Eourier mode k at 
each relaxation step n. 

Eigure[^ represents the relative update on the fields 


^(n+l) 


Jn) 


+ 


^(n+1) 


X(n) 


+ 


^rel 


^rel 


- 1 


(26) 


averaged over the spatial domain [0, U] x [—U, +E], as a 
function of the relaxation step n. In the relaxation algo¬ 
rithm, the relative updates on the fields decrease expo¬ 
nentially with the number of iterations before settling to 
some plateau. The algorithm described above therefore 
converges toward the solution of Eqs.( T2][l4 ). 



Relaxation iteration 


FIG. 1: Relative update on the fields Eq.(26), averaged over 
the spatial domain, as a function of the relaxation iteration n 


N 

f{u,v) = '^fk{u) cos 

k=0 



(i 4 , v) G [0, U] X [—y, +y] We can also have a glimpse on the numerical accuracy 

of the solution by evaluating the right hand side of 
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FIG. 2: Right hand side of Eq.(15) (in logarithmic scale) at 
the end of the relaxation algorithm (case of a solenoid with 
L = 0.1 and Cl = 1) 


the off-diagonal Einstein equation, Eq.(15), with the 
fields obtained at the end of the relaxation algorithm. 
This is illustrated in Figure in logarithmic scale, the 
rms of the residual is about 3 x 10“^ for field values 
of order unity (C/ = 1). The numerical error is more 
important around the sources of current and are mostly 
dominated by rounding errors in the evaluation of the 
elliptic functions that compose the magnetic potential 
Unr of the current loop and the solenoid [21]. 


We can now present the space-times curved by the 
magnetic fields of the current loop and the solenoid. 


B. Numerical Results 


Figures!^ and [^present the metric functions p and A as 
well as the relativistic part of the magnetic potential Urei 
for the current loop (Figure]^ and solenoids of different 
lengths (Figure [^. These plots have been obtained from 
the numerical resolution of Eqs.( p!2p4 ) with the bound¬ 
ary conditions Eqs.([T^ using the numerical method pre¬ 
sented in the previous section. 

Metric functions exhibit a peak at the loop location 
{u = l^v = 0) in Figure and we can notice the simi¬ 
larity between the solutions of the current loops and the 
solenoid of smaller length (Figure]^ first row: L = 0.1). 
Space-time is mostly curved along the radial direction 
{v = 0) inside the loop {u < 1) and becomes quickly flat 
outside of the loop. 

When the solenoid length is increased (Figure from 
top to bottom), we can see how the metric potential 
wells deepens and widens along the central axis. With 
longer solenoids, the total electric current involved is in¬ 
creased and so does the space-time deformation at the 
origin of coordinates. The peak in A near the solenoid 
location is also smoothed when the solenoid length is in¬ 


creasing (Figure]^ central panels from top to bottom). 
We also have that space-time deformation occurs mostly 
inside the solenoid and space-time becomes quickly flat 
outside the solenoid. 

The total magnetic potential a = a^r + <^rei is smaller 
than in the classical case since we have Urei < 0 on aver¬ 
age (Figures!^ and 1^ right panels). In general relativity, 
the electric current produces both magnetic and gravi¬ 
tational fields, leaving less energy to the former than it 
does in electromagnetism on flat space-times. 

These metric potential wells will produce light deflexion 
as well as gravitational redshift which will be maximal for 
a light source located at the origin of coordinates. For an 
observer located at spatial infinity (where gu 1) and 
a source located at the origin of coordinates, the gravita¬ 
tional redshift is simply z = 1 — exp(—p(0,0)/2), and is 
of the order of magnitude of the magneto-gravitational 
coupling Cj (see also Fig. |^. The precision achieved 
by optical lattice clocks in the measurement of a tran¬ 
sition frequency is of the order 10“^^ [12]. Achieving 
such a gravitational redshift with single-layered solenoids 
would require C/ ^ 10“^^, i.e. for an electric current of 
IkA^ n = 100 it would require a solenoid length of about 
lO^^m. Gravitational redshift therefore does not seem to 
be appropriate to detect gravitational fields artificially 
generated by coils with current superconducting technol¬ 
ogy. As we shall see further, Michelson interferometers 
with Fabry-Perot cavities will be more appropriate to 
attempt such detection of artificially generated gravita¬ 
tional fields. 

Fig. [5] illustrates how the gravitational redshift 
^ evolves with the magneto-gravitational coupling Cj 
Eq.([l^. For C/ ^ 1, the redshift 2 : varies linearly with 
Cl while a space-time singularity appears at the center 
of coordinates, p(0, 0) ^ —00 {z 1) when Cj oc. 

We can now focus on the relativistic effect of the devi¬ 
ation of light by an electric current, through the induced 
deformation of space-time by magnetic fields. This will 
be treated of the next section. 


IV. THE BENDING OF LIGHT BY MAGNETIC 
FIELDS IN GENERAL RELATIVITY 


In this section, we establish the general pattern of light 
geodesics in strongly curved space-times through numer¬ 
ical integration techniques. The geodesic curves x^{s) of 
neutral particles are solutions of 

cPx^ dx^ dx^ 

ds‘^ ds ds 

where s is some affine parameter. Introducing the metric 
ansatz Eq. this gives the following set of ordinary 
differential equations: 
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FIG. 3: Metric components and relativistic part of the magnetic potential for the current loop for Ci — 1 (left panel: gtt — 
— exp(p), central panel: grr = gzz = exp(A), right panel: Urei) 


(ft dp dt dr dp dt dz 

ds‘^ ^ dr ds ds dz ds ds 

d^z (? p_^dp(dt\^ ldX/dz\^ _p_x^P f dXdzdr 1 dX /dr\^ 
ds^ ^ 2 ^ dz \ds J 2 dz \ds J 2 ^ dz \ds J dr ds ds 2 dz \ds J 

d?r ^_^dp f dt\^ 1 dX f dz\^ r f dp \ f dp\^ dX dz dr l9A/dr\^ 

^ ~2^[Ts) + 2 ^ [d^J + 

d‘^p (2 dr dp dpdpdz ^ 

ds‘^ \r dr J ds ds dz ds ds 


(27) 

(28) 

(29) 

(30) 


Eqs.(27) and (30) can be 


dt 

ds 

dp 

ds 


directly integrated to give 



(31) 

(32) 


where we chose one integration constant such that the 
affine parameter s can be identified with the coordinate 
time ct at spatial infinity (where p ^ 0) and where the 
constant C is related to the angular momentum of the 


neutral particle. 

For null geodecics, we have that the tangent vector is 
light-like all along the geodesic curves g^y = 0 so 

that 


fi2 

e~P - 


= 0- (33) 


Putting this constraint into the remaining two geodesic 
equations Eqs.(28) and (29) gives 
















7 




1.0001 

1 

of 0.9999 
0)^0.9998 
0.9997 
0.9996 
0.9995 

0.9994 

10 




FIG. 4: Metric components and relativistic part of the magnetic potential for solenoids of different lengths for Ci — 1 (left 
panel: gtt = = exp(p), central panel: grr = gzz = exp(A), right panel: ttrei ; first row: L = 0.1, second row: L = 1, 

third row: L = 10) 

















FIG. 5: Gravitational redshift of a light source located at the 
center of coordinates and a receiver at spatial inhnity, as a 
function of the magneto-gravitational coupling Ci 


ds^ oz 2 \oz oz J \ds J 2 \oz oz J \ds J 

ds^ \or r J 2 \or or J \ds J 2 \or or J 


dX dz dr ^ 
dr ds ds 

dr\^ dX dz dr 
ds J dz ds ds 


(34) 

(35) 


In 

that 


the following, we restrict ourselves to planar trajectories in the (r, ^)—plane by setting (p = cst and (7 = 0 (such 
^ = 0). If we now set s = SI (and r = ul, z = vL) we finally obtain the following set of dimensionless ODEs: 


d^v ^ f dp dX\ f dv\^ f dp dX\ f du\^ dX dv du 

d?u LS f dp dX\ f dv\^ 1 f dX dp\ f du\^ dX dv du 


Since the metric fields and their derivati ves are all func¬ 
tions of {u^v) coordinates, solving Eqs.( 36|37 ) requires 
integrating both ODEs while carefully interpolating the 
numerical fields X^du,zX, p^du,zp at the current location 
point in each integration step. Eq.([3^ will be used as a 
constraint to validate the numerical integration. 

Fig. presents the deflected trajectories of a bun¬ 
dle of light rays incoming parallel from spatial infinity 
V + 00 . These trajectories have been obtained from 
the numerical resolution of geodesic equations in strongly 
curved space-times around loop and solenoids with ex¬ 
tremely large magneto-gravitational coupling Cj = 10 (or 
Cj = 1) so that the way light is deflected can be easily 
shown. The constraint Eq.(33) is below 10“^ along the 
trajectories shown in Eig. [6 Geodesics passing far-away 
from the current loop of the solenoid follow hyperbola as 


if photons were attracted by a point mass. However, null 
geodesics with close encounters can exhibit more sophis¬ 
ticated shape while passing through the magnetic device. 
In this case of strong coupling {Cj = 10 or (3/ = 1), the 
deviation of the light beam is so strong that magnification 
appears at some locations while large regions of the (i4, v) 
plane have been cleared of any light. These results can 
be applied to the gravitational lensing of cosmic string 
loops in the strong field regime. 


V. APPLICATION: GENERATION AND 
DETECTION OF ARTIFICIALLY GENERATED 
GRAVITATIONAL FIELDS 

We now investigate how far such deflexion of light 
could be detected with the present technology of super- 
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FIG. 6: Trajectories of parallel light rays coming from infinity (top left panel: current loop (Cj = 10), top right panel: solenoid 
with L = 0.5/ {Cl = 10), lower left panel: solenoid with L = I {Ci = 10), lower right panel: solenoid with L = 10/ {Ci = 1). 
Blue dots indicate the location of the current loop and the solenoids 


conducting electromagnets and high precision light-wave 
interferometers. The basic idea is to make interfering 
two light beams among which one has travelled in the 
space-time curved by powered superconducting solenoids 
and the other not. The shorter distances travelled by 
the light beam inside the powered solenoids will gener¬ 
ate a path difference between both light beams that will 
impact on their interference pattern as a result of a gravi¬ 
tationally generated phase shift. However, since the large 
electric currents that can be achieved with current super¬ 
conducting cables, roughly of order lO^H, will generate 
extremely weak space-time curvature, it will be necessary 
to amplify the signal by forcing light to perform numer¬ 
ous round trips in the artificially generated gravitational 
field. 

We can first write down the gravitational field equa¬ 
tions in the weak field limit. If we assume a^ei ^ ct^r and 
p, A ^ 1, we get that Eqs.( T2][T^ now reduce to 


{(duanrf - J (39) 


where + Ts reasons to be explained 

below, we will consider that the source of the magnetic 


fields is given by a set of stacked anti-Helmholtz coils 
(or multi-layered coils). The i—th anti-Helmholtz coil is 
constituted by two solenoids of radius /^ (/^ < /), length 
L, spaced by a distance D and carrying steady electric 
current of opposite directions. It will be necessary to pile 
up these (anti-)Helmholtz coils to produce a detectable 
space-time curvature by means of electric currents of or¬ 
der 10^A. The corresponding expression for the total 
magnetic potential Unr of n stacked anti-Helmholtz coils 
of radius k can be obtained from Eq.(ll) by 


n 

a„ = ^ (a-'(r, ^ + D/2- k) - a-'(r, ^ - D/2-, k)) 

(40) 

Corresponding boundary conditions can be obtained 
from those given in section H by linearly superposing 
the boundary conditions of single solenoids as allowed 
by weak gravitational field limit. The problem can be 
solved numerically with a standard spectral method (for 
instance based on Eourier decomposition, see section HI). 

Our proposed experimental set-up, shown in Eig. re¬ 
calls those of ground-based interferometers used to detect 
gravitational waves: it consists of a Michelson interfer¬ 
ometer whose arms are constituted by Eabry-Perot cavi¬ 
ties. One of these arms goes through a multi-layered anti- 
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FIG. 7: Schematic view of the proposed experimental set-up 

Helmholtz coil consisting of two stacks of superconduct¬ 
ing solenoids. As long as the electric current is switched 
on in the device, this one curves space-time and deflects 
light. Since space-time is slightly shrunk inside the coil, 
light trapped inside the powered coil accumulates phase 
shift as the round trips succeed each other. 

Following m, we can write down the phase shift due 
to the weak gravitational field on a Minkowski path 
7 of light 

A^=-f (41) 

2 

where is the (unperturbed) constant 4-wave vector 

of the wave front. We focus on light following the axis 
14 = 0 of the coil, with wave vector components = 

with u the frequency of light. Therefore, 
the phase shift along the axis of the coil for one trip is 
given by: 

TT 

= {p{0,z) - X{0,z))dz (42) 

where A = ^ is the wavelength of the light beam of 
frequency u and C is the length of the interferometer 
arm. An anti-Helmholtz coil configuration is better than 
a soienoid for the production of gravitationai phase shift. 
Indeed, phase shift is directiy related to the difference 
p — X which obeys the following PDE (Eq.([38|) - Eq.([^): 

V7p-A) = (43) 

where is the radial component of the magnetic field. 
This quantity, although vanishing on the axis of symme¬ 
try, increases more rapidly inside a Helmholtz coil than 
in the interior of a long solenoid, justifying the above- 
mentionned choice. 

As a matter of comparison, the phase shift induced by a 
gravitational wave passing by a Michelson interferometer 
is given by HU: 

‘2,71 c 

= — |/l|Ttrip (44) 


where \h\ is the amplitude of the incoming gravitational 
wave and Ttnp = 2£/c is the single round trip travel 
time of the light beam inside the arm (to be amplified 
through the multiple reflexions induced by the Eabry- 
Perot cavities). Currently achievable threshold for de¬ 
tection is for \h\ ^ 10“^^, jC ^ lO^m and A ^ 10“^m 
yielding A4> « 10“^^rad per trip. 

Let us now give an estimation of this phase shift for 
realistic experimental conditions of the setup presented 
above. We particularize the setup as following. We con¬ 
sider a set of 10 stacked anti-Helmholtz coils, each consti¬ 
tuted by two superconducting solenoids of same length 
L = 2.5m carrying opposite steady electric current of 
20kA (which is similar to CMS-class magnets [15]) spaced 
by a distance of D = 2.5m. The external solenoids have 
a radius of I = 5m and the 10 solenoid shells are chosen 
equally spaced between r = Im and r = 5m. The length 
of the interferometer arm has been chosen to C = 50m. 



v(l) 

FIG. 8: (p(0, 1 ’) — A(0, i’))/Ci (straight line) and magnetic field 
Bz{0,v)/{poNI) (dashed line) along the axis of symmetry of 
the anti-Hemlholtz coil (poNI = 4.4T). Dots indicate the 
position of the solenoids along the z—axis 

Figure 1^ shows the profiles of Bz{0^v) and p{0^v) — 
A(0, i;) along the axis of symmetry = 0 of the solenoids. 
The anti-Helmholtz coils generates a curvature of space- 
time that reaches its maximum at mid distance from each 
of the solenoids {u = 0, v = 0) where the magnetic field 
vanishes. The magnetic field is maximal at the center of 
the solenoids (at 2 ) = ±T)/2), and its magnitude reaches 
about 20T for the parameters chosen above. The phase 
shift Eq.([4^ can be integrated numerically for the re¬ 
sults of Figure]^ and we find Acj) ^ —1.56 x 10“^^ (for 
A = 514nm) per round trip inside the interferometer. 
If the experiment can be conducted long enough, for a 
time Texp, this phase shift will be accumulated. For 
two hundred days of duration, the accumulated phase 
shift reaches Acj) ^ —1.08 x 10“^^. This value is of the 
same order of magnitude to that of a gravitational wave 
signal [22], and could be detected by current technologies 
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developed for ground-based gravitational waves observa¬ 
tories m- 

The example given above is purely indicative and aims 
to show that a detectable phase shift could be produced 
by present day technology. We can give a simple order 
of magnitude and lower bound of the accumulated phase 
shift produced by posing p — A ~ C/ on the axis of sym¬ 
metry so that the phase shift Eq. (42) gives 


« TTjCl X 

^ ^bounce 


TT —CT 

'--^exp 


(45) 


where Texp is the duration of the experiment and tbounce 
is the time taken by light to produce a bounce inside 
the interferometer so that Texp/^bounce = ^bounce is the 
number of bounces inside the Fabry-Perot cavity. 

To conclude this section, we emphasize that the gen¬ 
eration and detection of artificial gravitational fields by 
strong magnetic fields is within experimental reach but 
requires large multi-layered superconducting magnets 
powered during dozens of days as well as hundred meters 
long Michelson interferometers with Fabry-Perot cavi¬ 
ties that could achieve the same sensitivity than ground- 
based gravitational wave observatories but in presence of 
intense magnetic fields. As in gravitational wave obser¬ 
vatories, a long optical path is crucial for the detection 
of the very weak gravitational perturbations. However, 
the wavefront of an astrophysical gravitational wave is 
much longer than the perturbations of space-time gen¬ 
erated by electro-magnets a decameter large. Therefore, 
a kilometer-wide interferometer can fit into the gravita¬ 
tional wavefront coming from astrophysical sources and is 
necessary to capture the wave during the time it passes 
through Earth. Using a kilometer-large interferometer 
for the detection of the space-time curvature induced by 
superconducting coils would require kilometer-large mag¬ 
nets. This is not necessary since, at the opposite of the 
detection of gravitational waves, space-time deformation 
by electro-magnets is maintained as long as the magnetic 
field is present. The amplitude of this space-time defor¬ 
mation is extremely tiny, of order of C/, which requires 
to trap the light long enough inside the space-time defor¬ 
mation to accumulate enough phase shift for detection. 
Although experimentally challenging, such a detection 
would open the path to a new class of laboratory tests of 
general relativity and the equivalence principle. 


VI. CONCLUSIONS 

The generation of artificial gravitational fields with 
electric currents could be in principle detected through 


the induced change in space-time geometry that results 
in a purely classical deflexion of light by magnetic fields. 
This effect does not invoke any new physics, as it is a 
consequence of the equivalence principle. Although very 
weak, we have shown that this effect could be detectable 
by a twofold experimental setup. On one hand, it in¬ 
cludes stacked large superconducting Helmholtz coils for 
the generation of the artificial gravitational field. On the 
other hand, the detection would be achieved by highly 
sensitive Michelson interferometers whose arms contain 
Fabry-Perot cavities to store light into the generated 
gravitational field. In an appropriate experimental set¬ 
up, the amplitude of the phase shift accumulated dur¬ 
ing the bouncing of light in the curved space-time gener¬ 
ated by the magnetic field would reach in a few months 
the level of an astrophysical source of gravitational wave 
passing through ground-based GW observatories. 

We claim that such detection would open new eras in 
experimental gravity and laboratory tests of general rel¬ 
ativity and the equivalence principle. These tests, al¬ 
though concerning the weak field regime, will have the 
particularity of focusing exclusively on the coupling be¬ 
tween gravitation and electromagnetism. Future theoret¬ 
ical works should focus on extending the present study 
to alternative theories of gravity to explore how far they 
would depart from general relativity. 

Such a detection of the space-time curvature generated 
by a magnetic field in laboratory would constitute a ma¬ 
jor step in physics: the ability to produce, detect, and ul¬ 
timately control artificial gravitational fields. And would 
this technology be developed, it could lead to amazing 
applications like the controlled emission of gravitational 
waves with large alternative electric currents. Gravity 
would then cease to be the last of the four fundamental 
forces not under control by human beings. 
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